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The optical fiber drawing process is considered in its totality — from 
source to forming zone to draw-down region and take-up end — as a prob- 
lem in fluid dynamics. Fiber draioing of most glasses is dominated by 
viscous stresses, surface tension effects, and quenching rates. This con- 
trasts with the drawing of textile fibers, where other fluid properties and 
non-Newtonian effects can play important roles. Preliminary time-in- 
variant "base flow" models are developed for glass drawing, using the one- 
dimensional, small-slope approximation of extensional flow. First-order 
sensitivities of these base flows to changes in operating conditions are 
examined via a stability analysis. Two important instability mechanisms, 
denoted as the tensile and capillary modes of dynamic fiber response, are 
discussed. Several follow-on objectives arising from this study are described. 

I. INTRODUCTION 

Stringent tolerances set on optical fibers used in communication 
systems have generated a need for understanding the fluid dynamics 
of the fiber drawing process. The responses of this process to various 
disturbances, especially those resulting in perturbations of the fiber 
diameter, are of interest. 

For steady-state drawing, one seeks analytic models that interrelate 
the draw-down ratio, draw force, flow rate, and some characteristic 
temperature. Such models serve two purposes : they predict the effects 
of changes in the operating parameters and can therefore be used in 
the control of industrial drawing processes; they also provide a "base 
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state" for the theoretical study of dynamic responses. Such responses 
could be due to a variety of physical disturbances; for example, 
mechanical vibrations, thermal transients, ambient gas flow, and even 
acoustic noise. All of these disturbances may produce variations in the 
diameter of the finished fiber. 

In studying the perturbations of liquid fibers, our philosophy departs 
somewhat from the more traditional one of stability analysis in textile 
engineering. Textile fiber studies, in general, strive to avoid fiber 
rupture and gross distortions of the thread line. (See, for example, 
Refs. 1 through 9.) In cases where continuous drawing of textile fibers 
is impossible, much emphasis is placed on the prediction and control 
of thread length, i.e., the filament lengths attainable between spon- 
taneous ruptures. 

The drawing of optical fibers takes filament continuity for granted. 
However, the ultimate optical application is sensitive to small diame- 
tral perturbations, far from rupture. Moreover, optical fibers are 
usually drawn at higher viscosities and draw forces than textile fibers. 
Starting from a much higher melt temperature, the subsequent viscos- 
ity increase due to quenching of the glass is much more severe than in 
polymers. Also, in its molten state glass is more nearly Newtonian than 
most polymers. 

Our purpose, then, is to model the mechanisms by which perturba- 
tions arise in the glass-drawing process and are frozen into the finished 
fibers. Based on this understanding, we hope to control fiber dimensions 
within the tolerances imposed by optical considerations (e.g. Ref. 10). 

The fluid dynamics of fiber forming involve a source flow, also re- 
ferred to as the forming zone, and a draw-down region. The forming 
zone is usually characterized by a rapidly contracting flow issuing 
from the bottom of a preform, or pulled from an orifice at the bottom 
of a crucible (Fig. la). In drawing from preforms, we may distinguish 
between a very sharply contracting configuration (Fig. lb), commonly 
encountered with laser heating, and a more gradual contraction that 
results from furnace heating (Fig. lc). We include in Fig. Id the case 
of an overheated preform, where the forming zone consists essentially 
of a liquid drop from which the filament is drawn. The different 
forming-zone configurations shown in Fig. 1 can imply fundamental 
differences in the flow field, as far as the steady-state and potential 
instabilities are concerned. For example, there is some evidence 11 to 
suggest that the overheated preform, Fig. Id, is capable of self-sus- 
tained oscillations at critical draw speeds, whereas, at noncritical 
speeds, the liquid reservoir in the pendant drop tends to absorb per- 
turbations coming from the take-up end. 

The forming zone makes a continuous transition to the draw-down 
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DRAW-DOWN ZONE 



(b) 



(d) 



FINISHING ZONE 



Fig. 1 — Fiber-drawing methods. 

region, where an essential part of the fiber contraction occurs. Typical 
draw-down ratios, defined as terminal velocity /source velocity, in this 
region run between 10 and 100. Most fiber stability studies place major 
emphasis on the draw-down region, because of its physical importance 
and mathematical tractability. However, the role of incipient perturba- 
tions in the forming zone cannot be ignored. The essential need for 
modeling this part of the flow field, if only by numerical simulations, 
is obvious. For completeness, we also recognize that a very small 
amount of fiber deformation occurs beyond the draw-down region ; but 
this falls within the visco-elastic rather than the fluid-dynamic regime. 
Much of the empirical evidence in textile and glass-fiber drawing 1 
suggests that fluid dynamics in the draw-down region is governed by 
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the so-called Trouton viscosity, a measure of the ability of fluid 
filaments to sustain tensile stress as a function of elongational strain 
rate. Indeed, if fiber breakage occurs in the draw-down region, it 
exhibits the neck-down and cohesive separation reminiscent of tensile 
test specimens. Much of the early modeling by Ziabicki 12 - 13 and sub- 
sequent efforts by Pearson, 614 Kas6, 6,7 and others are based on this 
notion, and for this reason we refer to these analyses as "tensile" 
models. 

In parallel, and possibly quite independent of the tensile mechanism, 
local flow conditions may exist where surface-tension-driven phe- 
nomena play a significant role. Such "capillary" responses could occur 
at the hot tip of the forming zone, where the viscosity is still quite low 
and the filament begins to take shape. Capillary models of filament 
response go back to Rayleigh's classical work, 15 ' 18 with subsequent 
extensions and experimental corroboration by, for example, Taylor, 9 
Tomotika, 17 ' 18 and Weber. 19 

The present paper serves several purposes: (i) it generates pre- 
liminary base-flow models for the draw-down region, using the one- 
dimensional, small-slope approximation of elongational flow; (u) it 
makes a preliminary assessment of sensitivities, through differences in 
base flow, to changes in operating conditions, such as take-up speed 
and quenching profile ; (Hi) it presents the tensile and capillary models 
of dynamic fiber response as fundamentally distinct mechanisms; it 
explores their applicability to different parts of the draw process by 
suitable modifications and extensions of existing theories; (iv) it pro- 
jects several follow-on efforts aimed at unified models of the steady 
and perturbed drawing process, viz. more realistic base-flow models, 
including heat transfer and two-dimensionality at the start of the 
draw-down region, and transient-response models, which account for 
these refinements in the base flow together with possible interactions 
between capillary and tensile mechanisms. 

II. REVIEW OF EARLIER WORK 

In this section, we discuss some of the literature on tensile and 
capillary stability analyses of liquid filaments. Table I relates several 
key publications and identifies their underlying assumptions and 
physical models. 

The first comprehensive study of tensile fiber models was undertaken 
by Ziabicki et al. 1-4 ' 12 - 13 ' 20 Motivated by the textile engineers' interest 
in potential instabilities and fracture mechanisms, the authors reviewed 
existing phenomenological evidence on filament "spinning."* They 



* This traditional terminology, which suggests twisting a fibrous material into 
strands, will be avoided henceforth as inappropriate to the drawing of liquid filaments. 
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Table 1 — Overview of tensile and capillary fiber models 



Tensile models, 
physical features 
discussed. 


Viscosity = 3 n, inertia, 
surface tension, 
quenching, finite 
draw-down. 


Viscosity = 3(i, 
stability analyses. 


Viscosity = 3 p, 

stability experiments. 


Pertinent literature 


Refs. 1-4, 12. 13, 20: 
analysis and experi- 
mental corroboration 
of relevant physical 
factors in base state 
and some dynamic 
perturbations. 


Refs. 5, 14 : dynamic 
analysis of attenuat- 
ing fibers, including 
the effects of viscos- 
ity perturbations and 
different base states. 


Refs. 6, 7, 21,22,25: 
experiments corrob- 
orating cyclic and 
transient responses of 
tensile fiber models. 


Capillary models, 
physical features 
investigated. 


Newtonian viscosity 
and surface tension. 


Inertia and surface 
tension. 


Newtonian viscosity, 
inertia, and surface 
tension. 


Pertinent literature 


Ref. 16 : cylindrical 
filaments neglecting 
inertia. 

Refs. 17, 18: adaptation 
to contracting flows 
and quenching effects. 


Refs. 16, 19: theory of 
inviscid capillary 
fluid cylinders. 


Refs. 16, 17, 18, 19: 
complete theories of 
capillary jets. Modifi- 
cation for quenching 
effects given here. 




Note : Numerous authors have conducted experiments over the years to dem- 
onstrate instabilities of capillary jets and filaments under isothermal, i.e., 
constant viscosity, conditions. 



also cite experimental evidence that when a filament is formed, say 
in the wake of a free-falling viscous drop, either a tensile ("cohesive") 
fracture or capillary separation may sever the flow. 

In Ref. 12, a suitable tensile theory is developed for the prediction of 
finite filament lengths, assuming the cross-sectional distribution of 
longitudinal velocities to be uniform and the steady-state flow to con- 
tract with small slope. Reference 20 accomplishes the same for capillary 
effects by an adaptation of Rayleigh's classical theory (as we have 
done independently for our own purposes). Ziabicki's subsequent 
publications include a thorough evaluation of experimental results, 
corroborating his predictions of filament lengths.* 

Ziabicki's comprehensive effort was followed by a series of papers 
by Pearson and others, aimed predominantly at tensile stability models. 
Starting from the simplest possible representation 14 — an isothermal 
filament under constant viscous tension — progressively more elaborate 
results were achieved by adding fiber quenching, inertia, surface 
tension, and gravity. A physical interpretation of this work, however, 
is difficult since the explicit features of diametral perturbation profiles 
along the fiber, their time dependence, and their sensitivity to proper- 
ties of the base state seem poorly understood. We will return to this in 
Section VII. Ziabicki's tensile flow models and the draw-resonances 



* Like other Western readers, the author has been somewhat late in fully recogniz- 
ing the significance of Ziabicki's work, much of which was initially recorded in Polish 
journals. Note, however, Ref. 4 for a more recent, comprehensive account. 
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predicted by Pearson et al. were corroborated by several experimen- 
ters. We cite the work of Kas6 21,22 and Donnelly and Weinberger 26 as 
examples in this area. 

As mentioned earlier, the study of surface-tension-driven perturba- 
tions contributes another, as yet separate, view of fiber stability which 
goes back to Rayleigh's classical work as presented in Refs. 15 and 16. 
Rayleigh himself studied several simplified cases : (i) constant viscosity 
plus surface tension (no inertia), (it) inertia plus surface tension (no 
viscosity), and (in) constant viscosity plus inertia plus surface tension. 
Each analysis yields an exact solution of the hydrodynamic perturba- 
tion equations for an incompressible liquid contained by a cylindrical 
boundary with surface tension. Simplifying assumptions are made only 
in solving the characteristic equations. 

Weber 19 showed that exact solutions of the characteristic equation 
differ little from Rayleigh's approximation. He also considered the 
effect of aerodynamic drag on the perturbed filament. Tomotika 17 - 18 
extended Rayleigh's model to allow for an ambient viscous fluid that 
surrounds the filament and is subjected to a steady elongational flow. 

For our own purposes, we need to modify Rayleigh's and Tomotika's 
work to reflect not only contraction in the base flow but also the viscos- 
ity buildup due to quenching. These are essential features of such 
"capillary" models of fiber drawing and are therefore listed explicitly 
in Table I. A display of diametral response profiles along the fiber and 
their dependence on wavelength and base-flow properties is given for 
comparison with tensile stability models. 

III. FUNDAMENTAL EQUATIONS AND THE ASSUMPTION OF 
ONE-DIMENSIONAL FLOW 

Let 2 = axial coordinate 

r = radial coordinate 

v = axial velocity component 

u = radial velocity component 

p = fluid density, assumed constant 

o- = surface tension, 

n = Newtonian viscosity, a function of temperature. 

In the cylindrical coordinate system (r, z), the Eulerian equations of 
mass and momentum conservation read : 

V. + U T + " = (1) 

P [Ut + UUr + WW.J = — + - 7V + ~j£ - - TB (2) 
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r i ii ^Ti dr rz 1 ,.,.. 

plv t + UVr + vd,1 = — + -T- + - r„, (3) 



where r r , r„ r„ denote stress components. In all other instances, the 
subscripts r, z, t denote partial derivatives. 

The constitutive relations lor an incompressible Newtonian fluid are 

,„ dv 

r r = - V + 2n- 

(4) 

1 o u 
T0 = —p + 2fi- 

r 



/dv . du\ 



where p is the pressure. It is one of the dependent variables, along 
with u, v, and the free surface configuration r = a(z, t). Equations (1) 
to (3) have the following boundary conditions : 

At z = : 

v(r, 0, = v (r, t) 

u(r, 0, t) = u (r, ( 5 ) 

a(0, = a (t). 

At z = L, the take-up position : 

v(r, L, t) = v L (r, t). (6) 

At r = a(z, t) : the kinematic condition 

da , da /n\ 

v = u— + — (*) 

dz dt 

together with tangential and normal surface-stress conditions, which 
we do not reproduce in detail at this point. (See Appendix A.) Once a 
solution of this boundary value problem has been found, the draw 
force at any cross-section follows from 



P(z, t) = 2ox<r/[l + (da/dz) 2 y + 2t / rudr. 

Jo 



(8) 



In particular, &P(t) = P(L, t) - P(0, t) and, to the extent that the 
solution for r t contains m, p, <*, the expression for AP depends on these 
fluid properties. 

The complete set of governing equations for the fiber drawing 
process includes an energy equation from which the temperature dis- 
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tribution T(r, z, t) in the draw-down region is determined. This, in 
turn, yields n as function of position and time. Strictly speaking, 
therefore, the heat-transfer equation is coupled to the fluid-dynamic 
equations. However, since we give it a separate, detailed treatment 
elsewhere, we uncouple it from this preliminary discussion and in- 
troduce /i(r, z), for nonisothermal draw-down, as a function presumed 
obtainable from some heat transfer model. 

The general axisymmetric, free-surface flow problem posed by (1) to 
(7) is a formidable one. However, for the purposes of an engineering 
analysis and to gain some basic insight, much headway can be made 
by taking advantage of the fact that \a,\ = 0(e) « 1 and \u/v\ 
= 0(e) « 1 throughout the draw-down region. This "small slope" 
assumption expresses the obvious fact that fluid flow in fiber drawing 
is essentially one-dimensional. The consequences of this kinematic 
feature are developed in Appendix A, taking advantage of the fact that 
we are dealing with low-Reynolds-number flow away from regions of 
strong relaxation in velocity profile. Specifically, we find : 

(i) v r = 0, i.e., "plug" flow throughout the draw-down region. 

(tt) T„ * € ( Tz) Tr , Tt)r/d 

(in) u a evr/a 

(iv) t t = re, uniform over the cross-section 
(v) T t = —a/a + 3juV ( , over the cross-section, 

where 3 n is the so-called "Trouton" viscosity for extensional 
flow in a liquid filament. 

These features of one-dimensional flow in the draw-down region sug- 
gest that mass and momentum conservation may be expressed con- 
veniently in terms of cross-sectional fluxes and stress integrals. Equa- 
tions of this kind may be obtained formally by integrating over the 
fiber cross-section: eq. (1) for volume conservation and eq. (3) for the 
axial momentum balance. Alternatively, we may derive these equations 
directly by taking a segment, of length dz, from the tapered axisym- 
metric filament as control volume. We obtain for volume conservation 

(o«i;).+ (a 2 )« = (9) 

and for axial momentum conservation 

p(aW) t + P (a 2 v) t - 3(a>*)* - <ra z = 0. (10) 

For an eventual comparison with the early work of Pearson et al., 
we also record (10) after the terms in p and a have been dropped. We 
have the simplified momentum equation 

o»m». = M(t), (11) 
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where M is an arbitrary time function. This equation merely states 
that the draw force is uniform along the fiber and varies only with t. 
For later reference, the time-independent (i.e., steady-state) versions 
of (9) and (10) yield 

a?v = const. = Q (12) 

and 

C = Sa 2 tiV. - P Qv + aa, (13) 

where C is a constant of integration, viz. the cross-sectional draw force 
reduced by the momentum flux. 

IV. STEADY-STATE SOLUTIONS FOR THE DRAWING PROCESS 

In this section, we consider solutions of the time-independent 
equations (12) and (13). These constitute steady-state representa- 
tions of the drawing process which are of interest for two reasons. First, 
they model the steady drawing operation and yield some insight into 
its controlling parameters, i.e., the dependence of draw force and 
draw-down profile on the draw-down ratio, the viscosity profile, fluid 
inertia, and surface tension. Second, they provide reference states on 
which to build dynamic response models for fiber-stability studies. In 
this context, such solutions are often referred to as base flow models. 

It is well-known in fluid-dynamic stability theory that detailed 
features of the underlying base flow can be quite important to the 
predicted dynamic response. Hence, it is necessary that we examine 
several base-flow solutions of the drawing process for the physical 
features they represent. 

Starting from (13), one observes that the first term on the right-hand 
side represents the viscous stress effect, the second fluid inertia, and 
the third a contribution from surface tension. We assume the following 
fiber dimensions and fluid properties at the start of draw-down: 

v = 10 cm/s 
a = 10 -2 cm 

v a = 100/s (e.g., Av = 100 cm/s, over Az = 1 cm) 
p = 2.5 gm/cm 3 
a = 200 dyn/cm 
H = 100 poise for soda lime glass 
= 1000 poise for fused silica. 

Note that the temperature at the interface between forming and draw- 
down zone is very dependent on as yet unknown fluid-dynamic and 
heat-transfer conditions in the forming zone. Therefore, the assumed 
values for m are rather tenuous. 100 poise probably represents a mini- 
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mum for soda lime glass, as might be expected in fiber drawing from 

crucibles. 

Given these data, the order-of-magnitude relations between terms 

in (13) are: 

Viscosity Viscous stress Inertia Surface Tension 

H (poise; 3juaV (dyn) pQv (dyn) aa (dyn) 



Soda lime glass 100 3 2.5 X 10~ 2 2 

Fused silica 1000 30 2.5 X 10 -2 2 

Thus, inertia effects amount to barely 1 percent of the viscous term 
under the most favorable circumstances, while surface tension can be 
important when n is near its minim um. However, physical evidence 
suggests rapid quenching of glass fibers in the draw-down region. 
Hence, inertia and surface tension effects soon become negligible as the 
temperature drops by several hundred degrees over the first few milli- 
meters of the draw-down region, causing fi to grow by several orders 
of magnitude (see, for example, Ref. 23). 

Disregarding, for a moment, the quenching effect that actually 
occurs in the draw-down region, we briefly consider an isothermal base 
state for two reasons. First, it permits an understanding of secondary 
physical effects such as inertia and surface tension, without being 
obscured by viscosity changes. Second, in later dynamic response 
studies, the isothermal base state serves as a basis of comparison for 
the stabilizing effect of the quenching that does occur. 

Letting the fluid properties in (13) be independent of z, we consider 
flow conditions such that inertia and surface-tension effects can be 
viewed as perturbations in relation to the viscous stress. Eliminating 
the radius a from (12) and (13) and nondimensionalizing according to 

*-v/vo, $ = z/L, ( )' = ^A 

we find 

tf ' - D* = - We& + Re*\ (14) 

where Vo = v at z = 0, L = length of draw-down region, and D = LC/ 
3a%v%ii , the nondimensionalized equivalent of C in (13). The inertia 
and the surface-tension terms are characterized by a Reynolds number 
Re = v Lp/3/i() and a Weber number, We = aL/3aoVono. The elementary 
solution for (14) with We = Re = is 

*«» - e^, (15) 

where we have used the boundary conditions 

* = 1 at r = 

* = E at f = 1 

and E = v L /v is the so-called "draw-down ratio." Note that D = InE. 
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Corrections to this simple base-state solution for small We and Re 
can be found by perturbation methods. A first-order approximation 
for We « D and Re « D is found by substituting (15) into the r.h.s. 
of (14) : 

•« = exp [ur + gj «-«•"»< - 1) + £ ( * : '" £ - 1} ] • (16) 






Resubstituting this into the right side of (14) for a second iteration, 
terms such as exp[e fIn *], were approximated by power series prior to 
quadrature with respect to f. Then, 

InE We + id 

2 » 



¥m = e finJ 



'(2) 



1 - We3C 



r jee_ 

L 2lnE ' 



We 
InE' 



I + \Re) ,~\ 
InE ' * J 



, Dnn \2Re me &We + \Re) "11 Q ~ 

+ Ee3C [z^'7^' M ' — m — >r Jr (17) 

where the expression for 3C[a, 6, c, d, f] is recorded in Appendix B. 
Numerical results from (17) are best presented in terms of the non- 
dimensional radius a/a = X. According to (12) 

X = *-*. (18) 

This has been plotted in Figs. 2 and 3 for a range of values in We and 
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Fj g 2— Perturbation solutions for base-flow profiles with increasing surface tension ; 
E = 10, Re = 0. 
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Fig. 3 — Perturbation solutions for base-flow profiles with increasing fluid inertia: 
E = 10, We m 0. 

Re. These results show a tendency for filament contraction to be de- 
layed with increasing surface tension and fluid inertia, a familiar phe- 
nomenon from more accurate base-flow models obtained by numerical 
methods. 24 

Now let us return to the more realistic case of variable viscosity. We 
denote it /to*?(f), where the second factor represents only the dimension- 
less dependence on f, with 77 (0) = 1. Then (14) becomes 



77*' - D* = - We& + Re* 2 . 



(19) 



Using an inverse approach, we can, for example, assume Re ^ 0, 
We = D = (where D = InE is no longer true) and * = e {lnE , as in 
(15). This yields 

Re 



v = r-s e * lnE 
InE 



(20) 



and constitutes an inertia-dominated base flow. Experimental data 
suggest that something like an exponential viscosity buildup along the 
draw-down region is a fair representation of quenching effects. Note 
that for this base-flow model, 

X = e-dnE/2)^ (2i) 

which has the disadvantage that X — ► for f ^> 1. To provide a finite 
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asymptotic value for X, we may assume a draw-down profile of the 
form X = f<r*t + g with / + g = 1, < /, g < 1, and solve for ¥ 
from (18) and rj from (19). 

We summarize the three elementary base-flow models obtained so 
far with the notation lnE/2 = a : 

Unquenched base flow ; Re = We = : 

X = e-«f, ¥ = e 2 <^, |»1. (22) 

Quenched base flow; fle ^ 0; We = 0; 



X — » for T -> co : 

X = e"°f, * = e 2af , 1? = &'*• 



(23) 



Quenched base flow; Re ?* 0; We = 0; 

X finite f or f — * °o : 

X = (/e— f + fl), * = l/(/e-r + <7) 2 , 17 = e a V(/e- af + 0>- (24) 

In addition to the inertia-dominated base flows (23) and (24), we can 
show that 

x = f<r<* + g, ¥ - ~(l " 2 f e ~ af ) . 1 " - g e< * W 

is a solution of the approximate fluid-dynamic and heat-transfer equa- 
tions, for small a, if we let Re = We = and assume 

M = Mfl efltt-r/ft), 

where T is the initial temperature and » 1. This represents a 
quenched base flow that is not inertia-dominated, in keeping with some 
of the perturbation equations discussed later on. Note that for each of 
these base-flow models da/dz = 0{aa /L). Since a = 0(1) but a /L « 1, 
this means that | da/dz | « 1 and confirms the basic assumption pro- 
viding for one-dimensional flow, as discussed in Section III. Note also 
that the viscosity profiles in models (23), (24), and (25), which reflect 
a cooling process along the fiber, are connected with the draw-down 
profile through the parameter a. This parameter is indicative of the 
quenching rate in 77, and also controls x'i the slope of the draw-down 
profile. 

The draw-force follows from any of these solutions by the obvious 
relation, in dimensional form, 

P = ira?T z + 2ax<x = 3iraV« + *&*< ( 26 ) 

If we neglect a and substitute one of the base flows, we find that P 
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depends on ao, Vo, no as well as the draw-down ratio and the quench- 
ing rate a ; all of which could be expected on physical grounds. 

Clearly, we could refine upon the viscosity profiles to be used in 
simple base-flow models. Since some experimental temperature profiles 
tend to show an exponential decay and the viscosity temperature 
relation for many glasses is of the form rj ~ e klT , we might consider 
functions of the form 

H — r + 8 exp (e 2af ) . 

However, such elaborations result in a loss of mathematical simplicity 
and usually lead to equations for ^ that require numerical integration. 
At that point, it seems more appropriate to solve the coupled heat- 
transfer/flow problem by numerical means. This has been done and is 
documented elsewhere. 24 

V. PERTURBATION. EQUATIONS FOR TENSILE STABILITY MODELS 

In this section, we develop the first-order perturbation equations 
necessary for a linear stability analysis of tensile fiber models. Let the 
first-order solution of (9), (10), and (11) be denoted 

a = o(2)[l + d(z, t)2 

v = t>( 2 )[l +0(M)] (27) 

- v = K*)[i + '(M)J\ 

where a (z), v(z), and v(z) represent radius, velocity, and kinematic 
viscosity for a suitable base state, in dimensional form for the time 
being, and d, v, P are dimensionless first-order perturbations of these 
quantities. 

Substituting (27) into (9), the first-order variation of the continuity 
equation reads 

*. + -*« + !*. = 0. (28) 

Similarly, (10) without the surface-tension term becomes 
. , . /*.,«. 2v\ 1 . , A /2v z 2v\ 

= -^-Kt)> (29) 

and from (11), the momentum equation without inertia terms, 

2d + , + + 1$, = ^), (30) 
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where we have used tfw, = Po = const, from (11) and d(t) is a per- 
turbation of M(t). Since the examples in later sections use the base 
states of Section IV, it is convenient to nondimensionalize the space 
and time variables as 

a! = I "-£ - T (31) 

and use the notation 

*Li _ ( )' and *L± - ( )'. (32) 

db dr 

The specific versions of (28), (29), and (30) now develop as follows. 
With the base states (22) or (23), eq. (28) becomes 

d' + W + de~ 2i = 0- ( 33 ) 

Neglecting inertia in the momentum equation, (22) and (30), for the 
unquenched inertialess base state, yields 

2d + H + P + W - <p(t), (34) 

where <p is an arbitrary function of r. On the other hand, the momen- 
tum equation (29), using the quenched inertia-dominated base state 

(23) leads to 

0" _ 4p _ 4dV-2» _ 2ije~ 2i = -2P' - 40. (35) 

Turning now to the inertialess, quenched base state (25), we revert to 
(28) and (30) as basic equations, but reserve the substitution of v and 
v t from (25) for a later time. 

For some of the examples treated in later sections, it is convenient 
to eliminate v from (34) and (35) by means of (33) and similarly from 
(28) and (30) for base state (25). The resulting equations for d are 
recorded in Appendix C for later reference. 

Conversely, the boundary conditions for some problems demand an 
equation in v. This is the case with steady-state responses to changes 
in the takeup velocity, which we treat in the next section. If (27) is 
used in (12) and we let Q - Q(l + d), the first-order variation of that 

equation yields 

2d + v = 4. (36) 

Next, consider (13), where we neglect surface tension and perturba- 
tions in v. As noted before, the constant C may be interpreted as a 
force parameter, carried from the forming zone to the draw-down 
region. Taking the first variation of (13), letting C = C(l + 6), and 
eliminating d by means of (36), 

*-£—&'+«(£-?)■ (37) 



OPTICAL FIBER DRAWING 1025 



Nondimensionalizing the base-flow variables in this equation according 
to Section IV, we have 

Re* . k . ( Re* ¥ f \ . ,_ Q . 

where 

LCc 



k = 



SdoVom 



Equation (38) is designed to yield changes in the steady-state velocity 
profile, v, as a function of $ and k, which are perturbations of the volume 
flux and force parameter of the base flow. 

In summary, this section has developed first-order perturbation 
equations based on the continuity and momentum equations of 
Section III. The perturbation equations were given in terms of & or v 
as needed for the steady-state and dynamic-response studies to be 
pursued in Sections VI and VII. 

VI. PERTURBATIONS OF THE BASE FLOW 

In this section, we use the time-independent versions of perturbation 
equations derived in Section V to display changes in several base-flow 
solutions due to shifts in such steady-state parameters as the boundary 
values and the viscosity profile along the fiber. Since these parameters 
are often accessible to control in real fiber-drawing processes, their 
effects on the steady flow are of operational interest. Obviously such 
effects could be determined by differencing neighboring base-flow solu- 
tions in the control-parameters space ; however, exhibiting the changes 
(analytically) as first-order perturbations can yield useful insight for 
the design of feedback controls. 

We start by examining the response of steady state (24) to a change 
$l in take-up speed. Substituting the appropriate base-flow expressions 
for * and t\ into (38), one has 



where 



v' - Re e — v = ke~ h X + (Re - 2/) %- 4, (39) 



a = «f ( )'-^ ■- /e- a +ff 



and a factor of 1/a has been absorbed in Re and A;. The boundary con- 
dition for (39) is 

v = v L at 3 = d L = a. (40) 

We find 

2(l + /)^\i2e V K%1) 
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with 



V - h + {' + We) «] ** + 2(TT7) (2 + Be//) ' 



where X L = X(2x,). Given (41), the corresponding changes in fiber 
radius d follow from (36). 

Additional features of interest with (41) are v at z = and perturba- 
tions of the draw force at either end. Thus, 

*-•<»- r-j^ + dJ-i)*. (42) 

For the draw force, 



which yields 
From this, 
and 



P = P(l + p) = 3irna%, 
*-* + ?« -* + -2JT 

Po = q + ~2T + 2(1+7) 



(43) 



^M + ^x^ + ^TTT)' ^ 



the changes in draw force, as function of i>l, 4, and &. 

The solution (41) is of little more than conceptual value as long as $ 
and k are unknown. Recall from (38) that these parameters represent 
perturbations in the integration constants of the base-flow solution for 
the draw-down region : $, a change in the volume flux, and k, a change 
in the force parameter of the momentum equation. Such changes must, 
in general, be expected to enter from the forming zone when the steady 
state is altered due to iiz.. 

Fortunately, ^ = for drawing from a preform that is fed at a 
constant rate. However, the exit flow from a crucible (Fig. la) does 
not provide such a simple condition. If operating at a low head in the 
reservoir, we would expect the entrance flow into the orifice to be 
affected by changes in the take-up speed. In neither case does there 
exist an obvious condition for the force parameter k [i.e., c and C, see 
eq. (38)]. 

A theory of the forming zone should be able to relate vo and po on 
the one hand with $ and k on the other. Given such relations, these 
would combine with (42) and (43) to determine 4 and k in terms of 
vl, and hence v, a, p as functions of vl- Depending on the different 
situations depicted in Fig. 1, the relations of v , po vs $, k in the forming 
zone could vary considerably. In some cases, an understanding of the 
complex fluid-dynamic and heat-transfer processes of the forming 
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zone (see Fig. Id) would seem essential for a satisfactory representation 
of speed-diameter-force relations at the take-up end. 

In the remainder of this section, we examine the sensitivity of steady- 
state flow in the draw-down region to changes in the viscosity profile. 
Such changes can be viewed as consequences of perturbations in the 
heat-transfer mechanism. Our primary purpose is to determine re- 
sponse amplitudes for d, the perturbation in draw-down profile, as 
functions of i>o, an amplitude parameter of the viscosity perturbation. 
Once again we encounter the problem of assuming reasonable boundary 
conditions at 2 = without a dynamic model of the forming zone. 

We consider two different cases: the draw-down response without 
inertia effects imposed on the unquenched, inertialess base state (22), 
and the response with inertia effects imposed on the quenched, inertia- 
dominated base state (23). In particular, we shall be working with 
eqs. (76) and (77) after deletion of the time derivatives. 

In both cases, we consider a viscosity perturbation of the form 

P = Pode-y*, (45a) 

which represents a distribution of arbitrary amplitude and spread, 
determined by Pq and 7, respectively. A family of such functions is 
displayed in Fig. 4 for 1 ^ 7 ^ 6. Note that for a given value of 7 

r 2 - - 1 . ( 45b ) 

"o max ye' 

the peak viscosity perturbation, normalized w.r.t. P . 

Let us consider various boundary conditions that may be applicable 
to solutions of this problem. If we assume that the fiber is drawn from 
a preform with constant feed and take-up conditions, an obvious 
boundary condition is 

d(d L ) = 0. (46a) 

In view of $ = 0, it seems reasonable to assume that the forming zone 
will respond with p(0) = 0. t)'(0) can be eliminated from these two 
conditions to yield 

a'(0)-^d(0) = 0. (46b) 

If three boundary conditions are needed, we take 

d(0) = 0, (46c) 

and then, according to (46b), also &'(0) = 0. Finally, for a fourth con- 
dition, let 

d'(d L ) = 0, (46d) 

which implies that p(3z.) = 0, i.e., no perturbation of the draw force 
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Fig. 4 — Profiles of the viscosity perturbation, v/v = de~">*. 



occurs at the take-up end. Conditions (46) will be invoked as necessary 
with increasing order of the perturbation equations. 

We start by considering the inertialess perturbation eq. (77). We 
have 



a = - -° e-y* ( - + A + Ba + Bt. 



(47) 



Results from (47) have been plotted in Fig. 5 to represent the effect 
of viscosity perturbations on the unquenched base state. Note the non- 
monotonic evolution of these curves with increasing y. 

Next, we examine the corresponding results, including effects of fluid 
inertia. Integrating (76) after deletion of the time derivatives, the 
quenched inertia-dominated base state (23) leads to 



6!" + 2d" - 4d' -8d= v" + 4/J' + 4/5 + C. 

C is a constant of integration. 

Substituting (45) into (48), we obtain 



(48) 



dOO = 



-Pt 



(7 + 2) 



o~i i 



(•+7*l)~ 

+ B^ + {Bi + B 3 d)e~ 2i + B 4 . (49) 
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Fig. 5 — Draw-down response to the viscosity perturbations of Fig. 4, without 
fluid inertia, for unquenched base flow. 

Figure 6 shows (49), the effect of viscosity perturbations on the 
quenched, inertia-dominated base state, over the range 1 ^ y ^ 6. 
Note that the presence of quenching significantly alters the evolution 
of response curves in Fig. 6, which becomes monotonic with y. Typical 
amplitudes in this case are about one-half to one-third as large as for 
the unquenched base flow (Fig. 5). 

The main inference to be drawn from a comparison of Figs. 5 and 6 
is that typical peak amplitudes for d/P are reduced significantly due 
to fluid inertia and quenching. 

By way of specific example, we consider results for y = 3.5: 

Maximum for viscosity perturbation (Fig. 4) P/P = 0.110 

Maximum for response without inertia on 

unquenched base state (Fig. 5) d/P = 0.018 

Maximum for response with inertia on 

quenched base state (Fig. 6) d/P = 0.005 

The latter case, which represents the more realistic model, also predicts 
the lower response amplitudes. Thus, for example, a 10-percent de- 
parture from the nominal viscosity profile would cause only a 0.5- 
percent departure from the draw-down profile. 
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Fi g _ 6— Draw-down response to viscosity perturbations, with fluid inertia, for 
quenched base flow. 

In summary, the first-order results given in this section yield a 
qualitative indication of base flow responses to time-invariant changes 
of the boundary conditions and of the viscosity profile. As stated in 
Section IV, a more satisfactory treatment of viscosity effects may be 
achieved by numerical integration of the base-flow equations, which 
introduce heat-transfer perturbations through the energy equation. 
However, the question of realistic interface conditions between the 
drawn-down region and forming zone remains open until the latter is 
included in our model. 

VII. THE DYNAMIC RESPONSE OF TENSILE FIBER MODELS 

We turn now to the dynamic response of tensile fiber models for the 
unquenched and quenched base flows; i.e., we address solutions of the 
equations in Appendix C, including the time-dependent terms. Un- 
fortunately, the formulation allowing for fluid inertia, eq. (76), does 
not lend itself to a simple solution. We therefore seek what preliminary 
insight can be gained from solutions obtainable with (77) and (78), 
i.e., by neglecting inertia in the perturbation equations. 

Starting with (77), which represents perturbations of the un- 
quenched base flow, the operator on the left-hand-side suggests a gen- 
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Fig. 7 — Tensile fiber model : frequency response of surface perturbation at 3>l = 2 
for base state [eq. (22)]. 



eral solution of traveling waveform 

&' = Sir + \e~™) 



and 



* - / fir + h<r n )dh + *(t), 



(50) 
(51) 



where $(t) is an arbitrary time function. This solution may be used 
to represent radial excitations at the source or take-up end or to satisfy 
boundary conditions in the presence of a particular solution. In the 
former case, where &(0, t) = sin cot, we reconstruct a solution by 
Pearson and Matovitch 14 of the form 



d(d, t) = Ai(3) sin cot + A i {'b) cos cot 
with the terminal response amplitude 

Ait*, n) = M(«#, t L ) + Alia, d L )y. 



(52) 



(53) 



This is normalized with respect to d , the amplitude of radial perturba- 
tions at 3 = 0, and plotted, for later comparison, in Fig. 7 as a function 
of co, for 3i = 2. It shows a series of response peaks presumably due to 
the absence of quenching from the base state (22), used in (77). 
These response peaks are commonly referred to in the literature as 
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"draw resonances." The author takes exception to this term since it 
implies the existence of natural frequencies contingent on the inter- 
action of system inertia and some restoring force, neither of which is 
obvious in the present model. Nevertheless, experimental evidence con- 
firms the occurrence of highly amplified responses near some of the 
"critical" frequencies predicted by Pearson's model. Typical radial 
perturbation profiles for d(3, r) are shown in Fig. 8 for w = 100, at 
T = 0, the start of a period, and T = 0.25, its quarter-point. They 
illustrate the spatial amplification of surface perturbations occurring 
along the draw path. 

Since the direct physical realization of radial perturbations at 
3 = may be difficult, we now examine the effect of viscosity perturba- 
tions that are convected along the fiber as a consequence of fluctuations 
in the heat source; i.e., 

Substitution into (77) yields the inhomogeneous equation 

(-?- + e»-f\d'= -g(r + h<r 2i ), ( 55 ) 

\dr dd ) 

where the dot is also used to designate differentiation with respect to 
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Fig. 8a— Surface perturbation for base state [eq. (22)] with w = 100, %l - 2, 
at T = 0. 



OPTICAL FIBER DRAWING 1033 




0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

DIMENSIONLESS DISTANCE g 

Fig. 8b— A at T = 0.25. 

the compound space-time argument. If we substitute 

6! = h(r + %e- 2i )e~ 2d , 

this leads to h = — $g ; hence, the particular solution 

& - -&(r + *e-'») = -*#. (56) 

This suggests that any space-time history of viscosity changes given 
in the form (54) translates into changes of fiber radius by the factor 
— \. It is a traveling-wave type of response only in the sense that it is 
convected with the moving fluid.* 

To take a specific example, consider a periodic viscosity perturba- 
tion, as would be caused by misalignment of the rotating laser beam 
used in heating the preform, 

= sin co (r + 17), (57) 



* The negative sign may appear surprising at first. One notes, however, from the 
continuity equation (33) written as 

(e~ 2d d/dT+d/dd)d = -\v', 

that & of the form g(r + je -28 ) leads to b' = 0. This means that such a form of radial 
perturbation can travel with the flowing fiber without perturbing the local velocity. 
Further, if a constant draw force is to be maintained, the expression 

P/p = irtfw' 

shows that, with v' unperturbed, an increase in v requires a decrease in a of half this 
magnitude. 
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where rj = \e~ 2i } yielding the particular solution 

dp = — h sin cot cos m ~ h cos cot sin coi?. (58) 

The boundary conditions to be satisfied with the help of (51) are then 

»(0, t) = 0, v(d L , t) = 0, (59) 

<U0, t) = isinco(T + £). 

In this case, the perturbations are driven only by the variations in 
viscosity. 

An appropriate form for / is 

/(*, t) = /i cos co(t + v) ~ / 2 sin co(t + v) ( 6 °) 

so that 

<U&, r) = HAcoscoT - / 2 sincoT)[Ci]^ 2 

- K/i sin cot + U cos u>t)\_SC\% 2 + *(r), (61) 
where 

[cotf= f cos "' 

J Ml 



'ui; 



W 



dtt> 






and /i, / 2 are integration constants. After determining $(t), /i» and 
/ 2 to satisfy (59), once again a solution of the form (52) is obtained, 
where Ai and At are recorded in Appendix D. It is interesting to note 
that A (co, 2> L ) for this case, if normalized in terms of P and plotted as 
in Fig. 7 shows exactly the same response spectrum, but with half the 
amplitudes. Moreover, the profiles of radial perturbations along the 
fiber for this case strongly resemble the ones obtained for radial excita- 
tion at d = (Fig. 8). 

The sharp response peaks given by the above solutions at certain 
frequencies reflect the absence of quenching in the base flow (22). In 
search of some allowance for quenching effects, we consider two ad hoc 
modifications of the tensile fiber model: the base flow (25), which 
includes moderate quenching together with moderate draw-down, and 
base flow (23), which represents quenched, inertia-dominated, exponen- 
tial draw-down. 

Using eq. (78), which is based on (25), a solution of this perturbation 
equation is developed in Appendix D for the case of cyclic perturba- 
tions in the starting radius d(0, t) = sin cot. The resulting expression 
for d(f, t) is recorded in (88). It turns out that the term tri* sin cot 
exceeds all other contributions to the dynamic response by several 
orders of magnitude for all values of co and y of interest. Thus, the 
perturbation in the fiber radius is merely a shift in the exponential 
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draw-down profile of the base state, whose spatial variation is inde- 
pendent of oj and characterized by 7, the quenching parameter of the 
base-flow viscosity profile v = v e 7f . This perturbation oscillates in 
time with constant phase along the entire draw path. Typical profiles 
of & are shown in Fig. 9, again for T = 0, the start of a period, and 
T = 0.25, its quarter-point. A contribution from the nonexponential 
components in d(f, t) is only seen at T = 0, the zero-crossing of its 
sin cot term. These profiles show none of the spatial amplification evi- 
dent in Fig. 8, which confirms, at least qualitatively, the attenuating 
effect of the quenching process. Unfortunately, a limiting comparison 
between the two models is not possible as the quenching effect is made 
to vanish, since that also requires a vanishing of the draw-down in base 
state (25). 

If the quenched, inertia-dominated base flow (23) is employed in 
the inertialess perturbation equation (30), the resulting model is indeed 
subject to criticism as logically inconsistent. However, as a plausibility 
argument, we might suggest that inclusion of inertia in the base flow 
would at least give a qualitative indication of changes to be expected 
from a more complete allowance for inertia effects. The formal exercise, 
starting from (23) and (30), closely resembles the derivation of (79). 
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Fig. 9a — Surface perturbation for base state [eq. (25)] with g = 0.75, u = 20, 
y = 2, at T = 0. 
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Fig. 9b— d at T = 0.25. 

The resulting frequency response functions evolve quite clearly from 
curves such as Fig. 7 with considerable smoothing of peaks and valleys. 
This suppression of the response peaks, due to quenching, certainly 
agrees with experimental evidence. 

The main point to be made in this preliminary assessment of tensile 
fiber models is that their frequency response curves and surface per- 
turbation profiles bear little resemblance to the perturbations caused by 
surface tension, which we discuss in the next section. To generate 
more realistic response predictions for tensile fiber models, including 
inertia effects and heat transfer, we will have to resort to numerical 
means. 

VIII. THE DYNAMIC RESPONSE OF CAPILLARY FIBER MODELS 

We next inquire under what circumstances the well-known phe- 
nomena of surface-tension-driven perturbations on liquid filaments 
apply in the fiber-drawing problem. Indeed, there may be limited 
portions of the draw path, presumably near the hot tip of the forming 
zone, where the viscosity drops low enough for surface tension to be- 
come significant. At least for low-melting glasses, such as soda lime, 
this is a possibility, as born out by the comparison of essential terms in 
the base-flow equations of Section IV. We shall characterize this type 
of fluid-dynamic behavior as capillary fiber models. As we shall see, 
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their spatial response profiles under harmonic radial excitation are 
totally different from those of tensile fiber models. By all indications, 
these two kinds of filament response are fundamentally distinct phe- 
nomena. The ultimate question is to what extent they coexist and 
interact in a real fiber-drawing process. 

Rayleigh's classical theory of capillary jets is our point of depar- 
ture. 16 - 16 It shows that the growth rate of "varicose" perturbations 
(axisymmetric harmonic surface modulations) is given by 



where 



= ff (ft - 1) 

2na[k 2 + 1 - k 2 P (k)ri(k)2 ' 

m = real, the rate of growth 
a = fiber radius 

X = wavelength of the perturbation 
k = 2ira/\ 
I and I\ = modified Bessel functions. 

The denominator of (62) turns out to be negative for all k.* 

If this capillary response model is locally applied to a base-flow model 
such as (24), assuming that base-state parameters change negligibly 
over the wavelength X, we may use it to construct the dynamic re- 
sponse along a contracting fiber. Then the evolution of a small surface 
disturbance may be synthesized using the stepwise relation between 
displacement amplitudes at successive instants of time 



U 



"+1 



- Unf^ + mu), (63) 



where 



U n and U n+ i = peak amplitudes of sinusoidal surface perturbations 
at t n and t n +i, respectively 



= stepwise scale factor due to fiber draw-down 



On 

m = m at center of wavelength 
At = tn+i — t n . 



Our response simulation convects the end points of a given perturba- 
tive wavelength at their respective speeds, while computing local fiber 



* Note that (62) results from a simplification of the characteristic equation, 
neglecting inertia effects, which in turn precludes initial conditions on perturbative 
velocities. However, it can be shown that the quantitative effect of this approximation 
on m is trivial (Ref. 19). 
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Fig. 10 — Response histories for different values of a. 

properties from the base model (24). We use the notation: 

X = a/a * = v/v 77 = n/no 

U = U/U , ol = aa , X = Ao/ao, 

and 

2 = (s + g)/2a , 
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where s, q, are end points of the perturbative wavelength. The non- 
dimensionalized surface tension is assumed to vary as 

a = a^ at , 

allowing for possible changes along the fiber (small values of fi approxi- 
mate a linear dependence). 

Figure 10 shows a set of response histories in terms of log {_U/X} for 
g = 0.1, X = 100, a = 0.25, $ = 0, rjo = 300, and 0.0008 ^ a 
^ 0.0016, illustrating the build-up to different asymptotic levels as a 
function of d. 

Figure 11 shows typical profiles of surface perturbations along the 
fiber for X = 15.7, which corresponds to a frequency of 200 Hz. As 
expected, the varicose response consists of sinusoids whose wavelength 
is progressively stretched due to fiber draw-down and whose amplitudes 
are modulated according to a response history such as given in Fig. 10. 
Note that this behavior differs drastically from the tensile fiber re- 
sponse of Figs. 8 and 9, which are also driven by radial harmonic 
excitation at the origin. 

Figure 12 displays asymptotic response amplitudes for a range of 
initial wave numbers 2ir/X and several values of g, the terminal radius 
in units of a for the base flow. For fixed g and d, the response has a 
maximum in the neighborhood of X = 50. The low responses at short 
wavelengths are due to vanishing of the Rayleigh instability as X — » 2ir 
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Fig. 11 — Plot of deflected surface shapes for % = 15.7. 
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Fig. 12 — Amplification factor vs dimensionless wave number for a = 0.001. 

while for large X , hence large z , the perturbations encounter a rapid 
viscosity build-up along the fiber. Contrast these response curves with 
the ones for tensile models, and a fundamental difference is again 
apparent. 

The influence of the terminal radius g, is illustrated in Fig. 13. As 
expected with a in the denominator of (62), increased draw-down 
(decreasing g) will enhance terminal perturbations. The second curve 
on that figure indicates the effect of /3, the growth rate of surface 
tension. Finally Fig. 14 shows the decrease in response with initial 
viscosity, t?o and the increase with rising values of o- , the initial surface 
tension. 

Since the above simulation averages fiber properties over a perturba- 
tive wavelength and does not ensure continuity of perturbative surface 
velocities between time steps, it seemed appropriate to corroborate it 
by a slightly different model, due to Tomotika, 17 - 18 which is also germane 
to our situation. In Tomotika's study, filament contraction is effected 
by a surrounding medium subjected to extensional shear flow, as in 
some of Taylor's experiments. 9 This apparent difference in base flow 
and the need to let ambient viscosity approach zero for our purposes 
seems to limit the applicability of Tomotika's model to the fiber- 
drawing problem. However, as we shall see, it agrees quite well with 
our adaptation of Rayleigh's theory. 
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Fig. 13 — Amplification factor vs g and /5. 

Tomotika finds the following asymptotic expression for interface 
perturbations as a function of several model parameters : 



where 



log [tf/X]* = »|S| j%(k)dk, 



(64) 



£ = ratio of filament to ambient viscosity 
k = 2ira/\ = local wave number of the perturbation 
a = nondimensionalized surface tension 
\//(k) = a kernel that is detailed in the references. 
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The behavior of this expression is illustrated in Fig. 15 by plotting 
it against k for various values of £. These plots show that, in the 
absence of quenching, the stabilizing effect of ambient shear flow 
diminishes steadily with ambient viscosity while the maximum response 
shifts to higher frequencies. This behavior is altered significantly if 
exponential changes of filament viscosity and surface tension are in- 
troduced to represent quenching effects, similar to our modification of 
Rayleigh's analysis. Equation (64) then changes to 






(65) 




Fig. 14 — Amplification factor vs tjo and &q. 
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Fig. 15 — Amplification factor vs k for £ -» « (Tomotika's model with constant 
viscosity and surface tension) . 



where 



»-'~-«-#[(if-'] 



and <p(k) is another kernel detailed in the references. Plots of this 
expression in Fig. 16 show good agreement with Fig. 12 for large £ 
(and g = 0.1). Plots of (65) for £ = 10 6 and various values of the 
terminal base flow radius are given in Fig. 17. They show the same 
response of [_U/X~\ m for g — > as Fig. 13. 

In summary, these results indicate that the two capillary models of 
quenched fiber responses, obtained by modifying Rayleigh's and 
Tomotika's analyses, are essentially equivalent. Note again that none 
of the response curves, such as Fig. 17, bear any resemblance to those 
of tensile fiber models. 

An additional piece of insight into capillary response mechanisms 
comes from Weber's work. 19 He reproduces Rayleigh's analysis by a 
somewhat different approach and obtains an exact equation for m, as 
well as a simplified expression that agrees with (62). Weber shows that 
the small errors in (62) are essentially due to the neglect of radial 
components of the flow field. He demonstrates this conclusively by 
rederiving (62) from a one-dimensional representation (recorded in 
Appendix E) which captures all salient features of the capillary re- 
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(b) TENSILE "NECK-DOWN" 

Fig. 18 — Comparison of capillary and tensile instabilities. 

sponse mechanism. Thus, multidimensionality is not the criterion that 
distinguishes capillary from tensile response models. The fundamental 
difference seems rather to lie in the energy source from which the 
perturbations are fed: surface tension in one case and axial stress in 
the other. 

Figure 18 attempts to emphasize this distinction in a pictorial 
fashion: (i) In the capillary model, surface tension, by overcoming 
viscous stresses, tends to accumulate fluid from both directions into 
periodic "beads," ultimately pinching off individual droplets as the 
minimum energy configuration. A multiplicity of such separations can 
occur independently of each other, and the essential physical mecha- 
nism is equally as valid for stationary, uniform filaments as it is for 
contracting base flows, (it) In the tensile model, on the other hand, 
the "worst" among random localized constrictions, due to surface 
perturbations, causes a tensile stress concentration which further 
reduces the cross-sectional area and results in a single, "run-away" 
tensile separation. This is the familiar necking of any tensile test 
specimen. The tensile stress associated with draw-down in the base 
flow is an essential prerequisite for this mechanism. Surface tension 
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will indeed contribute to tensile resistance in the necked down section, 
but it does not fundamentally alter the separation process. 

Thus, the capillary mechanism is a contest between surface tension 
and viscosity, independent of filament draw-down, whereas in the 
tensile mechanism the viscous stresses, jointly with surface tension, 
attempt to resist the draw force. 

IX. SUMMARY AND CONCLUSIONS 

Several important observations result from the discussions in the 
preceding sections. 

(i) Optical fiber drawing differs from textile fiber "spinning" in 
several essential ways. The flow of glass in the forming zone 
and draw-down region is dominated by viscous stresses. 
Inertia and surface tension play secondary roles. (They be- 
come noticeable only for the lower-melting glasses and then 
only in limited portions of the flow field.) This contrasts with 
polymer fiber forming, where, in some instances, viscosities can 
be quite low and other effects may be of comparable importance. 
Moreover, glass fibers are quenched over a wider range of tem- 
peratures than polymer fibers. This, together with the extreme 
temperature-dependence of glass viscosities, causes viscosity 
profiles along the draw path to rise much more abruptly than in 
textile fibers. Finally, and perhaps most importantly, molten 
glass can be considered very nearly Newtonian, which is not 
true for most polymers. Given the above physical features and 
the small-slope assumption of gradual draw-down, we can 
justify one-dimensional base states as useful representations of 
steady flow in the draw-down region. 
(ii) Given a base-flow model, its sensitivity to changes in operating 
conditions, such as the take-up speed and temperature profile, 
can be estimated by first-order perturbations. It turns out that 
the draw-down profile is relatively insensitive to significant 
viscosity changes, assuming that interface conditions between 
the draw-down region and forming zone have been modeled 
correctly. For nontrivial forming zones, e.g., Fig. la, b, or d, it 
is difficult to make reasonable assumptions for these conditions. 
Since we lack a complete understanding of the forming zone, 
but expect its flow field to change with perturbations in the 
draw-down region, our results must be considered tentative. 
(Hi) Fundamental differences exist between the tensile and capillary 
models of dynamic fiber response. The tensile mechanism seems 
to prevail in most of the draw-down region. For low-melting 
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glasses, we conjecture that the capillary model may apply in 
the short and very hot transition between forming zone and 
draw-down region, where surface tension can sustain perturba- 
tions that subsequently propagate by the tensile mechanism. 
The interactions of these two phenomena and their relation to 
dynamic responses in the forming zone itself are presently not 
understood. In the following we amplify each of these points to 
some extent. 

The base-flow models we supplied in Section IV are admittedly 
qualitative. A heat-transfer analysis was circumvented in this pre- 
liminary study by assuming exponential viscosity profiles, suggested 
by qualitative experimental evidence. Exact solutions of the coupled 
one-dimensional momentum and heat-transfer equations are now being 
carried out to allow for different kinds of heating in the forming zone 
and various cooling mechanisms in the draw-down region. The resulting 
simulation will be able to provide more detailed operational trade-offs 
between steady drawing parameters. It will also assess the limited in- 
fluence of fluid inertia, surface tension, and gravity. Finally, this 
modeling effort presents a natural opportunity for experimental cor- 
roboration by suitably instrumented steady-state runs, using laser 
and/or furnace-heated preforms or crucibles. 

As an extension of one-dimensional base-flow models, radial-heat- 
transfer mechanisms should be simulated, leading to nonuniform cross- 
sectional viscosity distributions at the start of the draw-down region. 
These viscosity distributions must be input to a perturbation model of 
axisymmetric free surface flow which generates the nonplanar velocity 
profiles expected in the transition between forming zone and draw- 
down region. Ultimately, the detailed flow fields of forming zones such 
as Fig. lb and d may have to be simulated by discretization techniques. 
If properly combined, these efforts may, hopefully, result in a unified 
base-flow model that properly allows for interactions between the 
draw-down region and forming zone in representing steady-state re- 
sponses to changes in the control parameters of the draw process. 

Finally, as mentioned before, it appears that vastly different fre- 
quency response curves and longitudinal profiles of surface perturba- 
tions characterize tensile and capillary dynamic responses as fundamen- 
tally distinct physical mechanisms. (Note the intuitive distinctions 
given at the end of Section VIII.) They do not seem derivable, in proper 
relation, from some universal fiber stability analysis. The question is 
then, what must be done to develop them into parts of a realistic and 
unified dynamic response model. 

Since analytic solutions for tensile responses of nontrivial base flows 
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in the draw-down region are not possible, numerical solutions by modal 
analysis or space-time integration are being attempted. This is ex- 
pected to shed further light on the anti-intuitive response profiles ob- 
tained from the "inertialess" perturbation equations. Similarly, the 
capillary model applied so far to an elementary base flow may be im- 
plemented, by numerical means, for conditions representative of the 
transition between forming zone and draw-down region. Combining 
these extensions of the tensile and capillary response models, it may 
be possible to relate dynamic records of thermal or mechanical surface 
perturbations coming out of the forming zone to diameter variations in 
the finished fiber. 
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APPENDIX A 

Implications of the Small-Slope Approximation 

We briefly sketch the consequences of the small-slope approximation, 
| a, | « 1 and \u/v\ « 1, as they evolve from (1), (2), (3), and (4) by 
scaling arguments. Let 

v = Vo*, u = Uo<p, /* = HM, V = Pot, , QQ . 

a = a x, r = a £, z - Lf , and t = - t. 

Then the small-slope assumption amounts to 

K =0 /«o\ <<L 

| v \ LJ 

Substitution of (66) into (4) yields 

q^r pj£ t + 2 ^UMo» 

L |_ ntf) J 

/toWo [ PoL „ 1 A juot 



L T ' 

fioVo _ 

Tr 



(67) 



L I moVo J ^ 

where p £/Wo = 0(1) and will be omitted henceforth. Substituting 
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T " = ^T 



(66) and (67) into (3), 
pv L 



Mo L \voao/ J 



= (-i + Wth + ^| lh(*s + eV f )]f. 



= « 2 (-i + 2,* r ) f + i [&(¥« + eV f )]i. 



Considering the case pVoL/no = Re <K 1, 
Hence, 

(&!*«)« = 0. 

If no constraint is to be imposed on tj, we have 

¥ { = (68) 



and 

_ f 
' L 

Now, from (1), 



t T m = -j- ai<p t = 0[e(T„ r r , t»)]. 



*r + (£r)<*+i>/t>-0 

and because of (68) 

n = v/t = -i*r- (69) 

Then, 

fr = fe = _ t - „* f (7Qa) 

f» = - J nf^rr 

and 

f, 1 + 2if* f . (70b) 

i is determined from the normal stress condition at the fiber surface 
£ = X. In dimensional form, 

Substituting (66) and (70) into (71) and dropping terms of 0(e 2 ) and 
higher, as well as a„, one finds 

i = | - n*t (72) 

for all values of £, where a = aL/noVoao. With this result, f« from (70a) 
becomes 

#. - - | + 3^ f , (73) 

where 3»7 constitutes the "Trouton" viscosity. (The additional r^f 
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term, augmenting 2i^ f in (70), originated in fy of (71). This, in turn, 
is due to <ps of (67), the cross-sectional contraction that necessarily ac- 
companies the extensional flow of fiber drawing.) We also note from 
(72) and (70) that 

(l + 7,* r ) { = T rt - T flj = 0, (74) 

which agrees with the radial equilibrium equation if (66), (68), and 
(70) are used in (2) and we let Re <s£l. 
The tangential surface stress condition at £ = X reads 

a,(T r - t.) + (1 + a 2 z ) T „ = (75) 

and, if rendered dimensionless, has leading terms of 0(e) ; hence, it will 
be ignored. For completeness, we also note that the kinematic boundary 
condition (7), rendered dimensionless and time-invariant, yields 

* r _ _ 2Xf 

which is the time-invariant continuity equation (9). 

APPENDIX B 

Second-Order Perturbation Term for the Base-Flow Solution [Eg. {17)] 

0C[a, b, c, d, f ] = r (l + ~ + d + § + f ) + § (l + d + f ) 

+ Ee K 1 + d + ?) {ae " nE " 2be ~ (lnmn) 

+ (1 + d) (^-eW» E + 2abeV nE ' 2 " - £<r' lnE ) 

-(- ^L e 3tlnE _|_ ^2 g(3U£/2)f _ H_ g-(3Zr.£/2)f 

lo o ■'J 

+ c j (1 + d) \aF(lnE, f) + bF (- ^ , f)] 

+ £ F(2Zn£, r) + abF ( Z -^ , r) + ~ *•(-**#, f) 



T 2 



raGan£,f)+6G(-^,r)]}, 
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where 
and 



«*»-S[("?+5)--5] 



APPENDIX C 

Differential Equations for First-Order Dynamic Perturbations in the Radius, 6 

Using (33) to eliminate v from the first-order variations of the 
momentum equations, we obtain from (35), which reflects the 
quenched, inertia-dominated base state (23) with inertia effects in the 
perturbation equation, 

&"" + d-'"e- 2i + 2d"' - 4d-"e~ 2i - 2&"'(T a + 8d"e-** 

- 4d" - 8&' = W + 40" + P'". (76) 

From (34), the inertialess momentum equation, and the unquenched, 
inertialess base state (22), we obtain 

d" + d-'e- 23 = i>\ (77) 

Note in (76) and (77) that ( )' = d( )/dd and ( )■ = d( )/dr, 
where d = az/L and r = av t/L. Finally, from (28) and (30), the 
inertialess momentum equation with the quenched inertialess base 
state (25), we can find 

where ^ = v/vo as in (25), f = z/L and r = tvo/L. In this case, 

( y = d( )w. 

APPENDIX D 

Detailed Results tor Forced Dynamic Responses of the Tensile Fiber Model 

The detailed expressions for Ai and A?, resulting from (58) and (61) 
are 

Ai = | g [cos]^ + j A^gg X C(-Bi sin »/2 + 5 2 cos »/2) 

X [Ci^l. - (Bi cos o>/2 + £ 2 sin w/2) X [j5K]$ J 

Ai « | i [sin]^ 2 + 4e( ( ^ f "^ C ^ i) X [(Bi cos w/2 + B 2 sin «/2) 

X [CftlA + (-5: sin w/2 + 5 2 cos w/2) X [50$.]} , 
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where e = e 2d L and 

l«/2 



B 1 = [ Ci + ~ Si + cob]"' 

Si - £ Ci + sin 

Ze J M /2e 



(79) 



Recall that d = az/L. 

We now parallel the development of Section VII for the quenched 
base flow (25). Equation (78) with P = suggests a solution of the form 

d(r, r) = Re j/o jf expOU"' - T)] ex P r>(r - {)]#' 

+ 3> exp(l'coT — af)| , 

where 

$ « /'' #!- , /o = A + */», #o = *i + »*.. (80) 
Jo v\i ) 

Explicitly, 

d(f, t) = Ai sin ur + Ai cos ur, 
where 



Ai(J-) = /" f e a( f'- f) (/isina.^ - / 2 cos co£)cir' - $ 2 e- ar , 
Jo 

A 2 (D = /" e a(f '- f) (/icosco^ + /asinco^df + $ie-" f , 
Jo 



(81) 



describing the f -dependent phase and amplitude of the response. Note 
that f = z/L and t = tv /L, whereas in (79) b = af and r = a»rf/L. 
We use the same boundary conditions as in developing (52). For 
cyclic radial perturbations at f = 0, 

&(0, t) = sin cot, (82) 

together with the velocity conditions 

0(0, r) = 0(1, r) - 0. 

An expression for v is obtained by eliminating v' from (28) and (30) 

H= *d' + *d> -2d + 2G(r), (83) 

v v 

where G'(r) is an arbitrary time function and the differentiation symbols 
mean ( )' - d( )/dr, ( )' = d( )/d£. Substituting (81) into (83) 
we find 

j)(f, t) = Vi sin cot + V-i. cos cot 
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with 

-2o> ft 



Fl(f) = ^jT /o e " (r ~ r) (/' cos "* + /« sin «*)#' 

- 2 /" e-(f'-r) (/ x sin o>£ - / 2 cos «£)#' 

+ 2 ^[(* a + l) *,_£*,] 



2«J- 
+ -f (/x sin a>£ - / 2 cos «{) + 2Gi 



FiO") = I? / o f e a(f '- f) (/i sin »{ - /, cos »{)#' 

- 2 /" e«(f'-f) (/ x cos a>£ + /« sin w{)dr' 

-^■<|> *. + (£« + !)••] 



Ultimately, (82) leads to 



+ ^ (/i cos a>£ + /> sin «{) + 2G 2 . (84) 



where 



$i = 0, $ 2 = -1 

A = (fiiCi + B 2 C 2 )/(Bl + 51) / 2 = (/i5 2 - C«)/5, 

Gi = 1 + 77 (a + / 2 cos «£ — /1 sin a>£ ) (85) 

G 2 - — 77 (w + /1 cos u£ + A sin £ ), 
\ro 

£x = £ #(1) - K(X) + ft cos a»{i - 1 
B t .-«JC (1 )-ff(l)+ tiring, m 



*-(*« +1 K'- s- 1 - 



The subscripts 0, 1 denote evaluation at f = 0, 1 respectively, and the 
quantities H, K are denned as 



H(t) = /"%«<*'-» sin «&#' 
Jo 

#(D = /" e««"'-« cos «!#'• 
Jo 
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(87) 



With this notation, we ultimately get 

<*tt, r) = C/iffCT) - M'(f) + e- a f]sinajT 

+ C/i-K (f) + /«» (J")] cos tar. (88) 

APPENDIX E 

Weber's Derivation of the Capillary Stability Equation 

This appendix gives a simplified derivation of the stability equation 
for the capillary fiber model based on assumptions that are quite 
equivalent to those made for the one-dimensional tensile model. In 
fact, the rationale used here closely parallels that of Section III. 

The analysis proceeds in terms of equilibrium and continuity equa- 
tions, which we write in dimensional form for the entire filament cross- 
section. The perturbed surface radius becomes a + 5 and a given cross- 
sectional element is displaced by vdt along the fiber over the time 
increment dt. The radius of this element now becomes 

a + l + ^dt + fvdt. (89) 

dt aZ 

Since the last term is of higher order, it will be neglected. 

The constitutive relations are as in (4) and the derivation of an ex- 
pression for t, is quite similar to Appendix A. The main difference arises 
in the radial stress boundary condition, where, in distinction from the 
treatment of (72), the longitudinal curvature term cannot be neglected 
for varicose perturbations. Then the r.h.s. of (72) becomes 



(!-") 



(90) 
\a* / 

and, instead of (74), 

T , -*(£ + «") + 3m»'. (91) 

The continuity equation yields 

^ + t/ = 0. (92) 

a 

Now, combining (9) and (10), 

t' z = pi) -\r pw', (93) 

where the last term was apparently overlooked by Weber but seems to 
have little effect on the resulting stability equation. 

Substituting (91) and (92) into (93) and assuming surface perturba- 
tions of the form 8 = b*e mt cos kz/a, leads to the stability equation 

m> + 7n^4k* = 7 ^- 3 (l- k*)k\ (94) 

pa 1 Zpa 6 
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The unstable root, of interest here, is 
(3m)* 2 



m = — 



2pa 2 



+[^< i -^+^]' (95) 
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